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ABSTRACT 


Effects of diffraction and non-1 inear photographic emulsion characteristics 
on the performance of deferred electronic heteroayne moire def lectometry are 
investigated. The deferred oef lectometry is used for measurements of 
non-steaay phase objects where it is difficult to complete the analysis of the 
field in real time. The sensitivity, accuracy and resolution of the system 
are calculated and it is shown that they are weakly affected by diffraction 
and by non-linear recording. The features of the system are significantly 
improved compared with the conventional deferred intensity moire technique, 
and are comparable with the on-line heterodyne moire. The system was evaluated 
experimentally by deferred measurements of the refractive index gradients of a 
weak phase object consisting of a large KD*P crystal. This was done by photo- 
graphing the phase object through a Ronchi grating and analyzing the trans- 
parency with the electronic heterodyne readout system. The results are 
compared with the measurements performed on the same phase object with on-line 
heterodyne moire def lectometry and with heterodyne holographic interferometry 
methods. A good agreement was observed. Some practical considerations for 
system-improvement are discussed. 

In the second part of the project, an algorithm for reconstruction of 3-D 
phase objects from incomplete deflections data, has been developed. The 
algorithm is based on the inverse Cormack transformation. The advantages of 
the present reconstruction method are: 1. The measured deflections are readily 
used in the inversion transform; no numerical derivations have to be performed 
like in the case of interfercmetry. 2. ho iterations between the object domain 
and the projection domain are necessary. The algorithm has been evaluated by 
«a computer simulation; an asymmetric analytical density field around a 
circular opaque object was successfully reconstructed. 
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FIGURE CAPTIONS 

Fig. 1. Schematics of a conventional moire system. L.S. is a point light 

source (may be a focused laser beam). M is a parabolic mirror. G^ 
and G 2 are Ronchi rulings and S is a mat transparent screen. 

Fig. 2. Schematic of the setup for recording the phase object. <b is the 

refraction angle; R(x) is the distance from the phase object at which 
the refracted ray seems to emerge. 

Fig. 3. Schematic of the experimental setup for postanalyzing phase objects: 

G, Ronchi ruling; R.S. and T.S. , reference and test signals 
respectively; G.S., gating signal. 

Fig. 4. Intensity transmi ttance - exposure curves . 

a) For Agfa Gevaert 10E75 plate; the dashed curve is the third-order 
polynomial approximation, and b) For Kodak G49F plate; the dashed 
curve is the second-order polynomial approximation. 

E|_ and E Q are respectively the width of the linear portion of the 
curves and the exposure at which they start. 

Fig. 5. Calculated moire intensity vs. phase for various values of f. f = 0.0, 
0.02, 0.04, 0.10 , 0 .15, 0.23, 0.25, 0.33, 0.35, 0.40, 0.45, and 0.50 
for (a), (b), (c), (d), (e), (f), (g), (h), (i), (j), (k), and (1) 
respectively. In each figure, the upper two curves are tor 

0 

Agfa-Gevaert 10E75 plates prefogged with exposure E Q = 12 erg/cm 

? ? 

and with E 1 = 8 erg/cm and E' = 18 erg/cm for the solid line and 
dashed line, respectively. The lower curve is for Kodak 649F plates 
with E q = 0 and E' = 55 uJ/cm 2 . 
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Fig. 6 


Fig. 7 

Fig. 8 
Fig. 9 
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I 

I 


: x - derivative of refraction index for flow simulator. 

a) Deferred heterodyne moire 

b) On-line heterodyne moire. 

: x - derivative of refraction index for flow simulator - heterodyne 

holographic interferometry. 

: Object field and the geometry of the transform. 

: The object function used in the study of the algorithm. The 

functions is defined by Eq. 12. The segment r<0.6R Q is missing. 
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PART 1 PERFORMANCE OF MOIRE DEFLECTOMETR Y WITH DEFERRED 
ELECTRONIC HET ERODYf£ READOUT 


I. INTRODUCTION 

Moire def lectometry is a simple technique for optical testing of phase 
objects and specular surfaces, based on the moire effect-*- -4 . The method 
provides mapping of ray deflections caused either by a phase object or upon 
reflection from a surface. From this information the index of refraction 
field of the phase object or the quality of the reflecting surface are 
obtained. A summary of progress through 1984 is contained in the reviewing 
article by Kafri and Glatt^. 

A. Backgroun d 

A conventional moire def lectometer consists of two identical Ronchi 
gratings, a collimated light source and a diffusing screen attached to the 
output grating. A schematic description of the system is shown in Fig. 1. 
The gratings are shown separated by a distance a and their lines are rotated 
relative to each other by a small angle e. When a collimated light beam 
passes through the gratings a moire pattern is produced on the screen. The 
pattern consists of straight fringes perpendicular to the original grating 
lines, separated by a distance of p'*- 



where p is the pitch of the gratings. 

The moire effect for very small grating separation distance can be 

explained by pure geometrical optics. However, when grating separation 

increases, diffraction introduces undesired effects which limit the 

fS 7 

performance of the system ’ . To minimize these effects, the distance a 
between the gratings must be exactly one of the Fourier image planes of the 


- 2 - 


gratings namely 6 ’^ 



z =1 >2,3 ... . 


( 2 ) 


For these distances, the fringe contrast will be maximum. Wien 


A 


u 


+ I) £ 

2 ‘ x 


(3) 


the fringes will vanish completely. In Eqs. (2) and (3) x is the wavelength 
of the collimated light beam. 

It the collimated light is refracted by a phase object before it enters 
the deflectometer , the original straight parallel moire fringes will be 
distorted. d (r), the angle of refraction in the x-airection at a point 

A 

r = x,y, is related to the fringe shift 6h^ (r ) in the y-airection by 


"x (r > 


e sh y (r) P «h y (r) 
a = " p 'a 


(4) 


By measuring 6h^(r), d x (r) can be calculated. In turn, the x component of 
the index of refraction gradient may be calculated by using the equation* 


d x (r) » 


f z. 


an(r ) 

3 X 


dz 


( 5 ) 


or 


<M r) - J Zf d Z (6) 

x ' n^ z J 3x 

o 

where z q and z f are the boundaries of the phase object along the line of 
sight, n is the index of refraction which is related to the density p by 


n - I = 6 p 


(7) 
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where the coefficient 6 is the Gladstone- Dale constant, a function of the 

9 3 

wavelength x. For air, at a wavelength of x = 6328 A, G = 0.227 cm / gr. n^ is 
the index of refraction of the medium surrounding the object. 

B. Electronic Heterodyne Method 

Recently, the accuracy and sensitivity of the moire method, for on-line 

real time measurements, has been significantly improved by applying the 

8-1 0 

heterodyne technique for reading the deflections of the moire fringes 
The basic idea of the heterodyne technique is to measure the phase of the 
signal from a photodetector detecting the light transmitted through a 
travelling moire fringe pattern. The moving fringes may be accomplished by 
translating the output grating G„ relative to G. at a constant speed V in 
the x-direction. Khen the collimated light is not disturbed by the phase 
object (straight line fringes), the first harmonic of the photodetector signal 
is given by 9 

I (y,A*,t) = cos (nA* )cos[2n (-^ + + ^)] (8) 

if 

where a* = a/(p 2 /x) is the distance between the gratings in Fourier units and 
x/p and yo/p are phase shifts related to the relative translation and rotation 
of the Ronchi gratings. Here x is the x-directed offset between the lines of 
the two gratings, ft is the fundamental frequency 

ft = 2ir ^ (y) 

In deriving Eq . 8 the following assumptions were maoe: 

1. Paraxial approximation. 

2. Slow variations of the phase object over a pitch scale. 

3. Smal 1 e. 

4. Fresnel diffraction. 
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The effect of refraction by a phase object is to introduce a phase shift 9 
(r)A 

A<Mr ) = 2 t p (10) 

A procedure for measuring this phase shift is to record the electronic phase 
relative to that of a reference signal, with and without the phase object, and 
subtracting the two phases. Equations (4) and (6) are then used to compute 
the properties of the phase object. Since phase measurements can be 
accomplished for better than a degree, fringe interpolation accuracies better 
than 1/360 fringe are possible. 

C. Transient Phase Objects 

Mien unsteady phenomena occur, it is difficult to complete the analysis of 
the whole field in real time. In these cases, it is necessary to freeze the 
phase object and post-analyze it after the fact -1 ’ . The idea of applying 

the electronic heterodyne readout technique for the deferred moire fringes 
readout was recently demonstrated . In this case, the second Ronchi 
grating (see Fig. 2) is replaced by a photographic plate. The phase object is 
recorded on the photographic plate by exposing it to the collimated beam 
transmitted through the phase object and the first grating, thus the distorted 
image of the first grating due to refraction by the phase object is recorded. 
Readout is accomplished by placing the moving grating immediately in back of 
the processed photographic plate as shown in Fig. 3. When the combination is 
illuminated, a time varying signal is obtained and analyzed by the heterodyne 
technique in a manner similar to that described for the on-line case. 

The purpose of this work is to investigate, both analytically and 
experimentally, the effects of diffraction and non-linear response of the 
photographic film on the performance of the deferred heterodyne moire 
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def lectometry. The resolution, accuracy, and sensitivity of the method will 
be evaluated and compared witn those of the real-time moire deflection 
method. Some practical considerations for the most efficient use of this 
method will be discussed. 

In the next section, the theory of the method is presented and in Sec. Ill 
mapping the index of refraction of a weak phase object using this method is 
demonstrated. 

II. THEORY 

As mentioned above the measurements of the density field of a transient 
phase object involves three steps: Recording, moire fringe formation, and 
electronic heterodyne readout. First, the distorted shadow of the first 
grating is recorded on a photographic plate when placed in the position of the 
second grating - Fig. 2. The grating is illuminated with the collimated beam 
passing through the phase object. The photographic plate is later processed 
to a transparency, which later is referred to as the phase object 
transparency. The distortion of the grooves of G are mainly caused by three 
reasons: 

1. Refraction by the phase object. 

2. Diffraction by the Ronchi grid. 

3. Non-linear response of the photographic emulsion. 

This step is repeated when the phase object is absent, the transparency 
obtained is called the reference transparency. Even in this case the image of 
the grooves may be distorted as a result of diffraction and non-linear 
photographic effects. 

Next, the moire fringes are formed by placing the shadow transparency in 
front of the second Ronchi grating. When the combination is illuminated, 
moire fringes are observed. Since the shadow transparency and the Ronchi 
grating are in contact, the formation of moire fringes is a purely geometrical 
effect. 
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The readout step is accomplished by translating the grating (thus causing 
the moire fringes to move) and recording the phase of the signal output of the 
photodetector (Fig. 3) as described for the on-line case. 

In the course of calculating the output signal of the photoaetector we 
will follow the three steps described above. 


A. Intensity Transmission Distribution of the Transparency 

Assuming an infinite Ronchi grating with square wave line profile 
illuminated by an infinite beam, and assuming Fresnel diffraction, the 
intensity distribution of the beam passing through the phase object and the 
Ronchi grating at a distance a (where the photographic plate is located) is 
given by^: 


I (u , a* , 



+ 


- I ? 1n [ ^ 1)Uq cos[,(2n+l) 2 u t «] 
11 n=0 ^ 2n 


+ 


+ 4__ £ sin 2 [2Tt(2n+l )ug] + 8 

* 2 n=0 (2n+l) 2 n 


V s in[2*(2n+l )u£] 
n£m (2n+T] 


x ' Sin ^ 2 (2m^ I7 )U ^ cos ^hua*[ ( 2n+l ) 2 - (2m+l) 2 ])] (11) 

where 5 = x/p, u = and R(x) and d x (r) are defined in Fig. 1. In 

deriving Eq . (11) the paraxial approximation was assumed, which in this case 
means that k (x ) >> a,x. According to this assumption it is clear that the 
intensity in Eq . (11) is practically independent on the origin of 
coordinates. Note that when no phase object exists, u = 1 . The transmittance 
of the photographic plate is proportional to this spatially varying intensity 

p 

Generally it is assumed that there is a linear relation between IT! , 
the intensity transmittance of the photographic plate, and the exposure E. 

O 

However, the actual IT ( -E curves for photographic materials are always 
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p 

nonlinear to some degree. Figures 4a and 4b display the |TI -E curves for a 
Agfa-Gevaert 10E75 plate and for a 649F Kodak plate, respectively. In the 
figure, E Q denotes the exposure at which the linear region starts and E^ 
denotes the width of the linear portion of the curve. 

If the photographic plate is uniformly prefogged by exposure E Q , and 
illuminated by a pattern of amplitude E' (the total exposure is E Q + E‘ ) 

O 

then a moderate portion of the ITI -E curve can be approximated by a cubic 
polynomial , 

11,2 = R 0 + r i( E q + E ') + r 2 ( E o + E') 2 + R 3 (E q + E') 3 (12) 

Figure 4a illustrates the polynomial approximation for the Agfa-Gevaert plate 

with E q = 12 erg/cm 3 , 0 < E' <18 erg/cm 3 , R Q = 1.1847, R^ = -7.4385 x 10 3 
-3 

R 2 = 1.5 550 x 10 ', and R^ = -1.0726 x 10 In this exposure range, the third order 

term is small and may be neglected. The cashed curve in Fig. 4b is a quadrature 

2 

polynomial approximation for the Kodak plate; here E q = 0, 0 < E' < 53 y«3/cm , 

R q = 8.6668 x ID" 1 , R 1 = -5.3625 x 10 -3 , and R 2 = - 7.3440 x 10 -5 . Replacing 
E' in Eq. 12 by I, given in Eq. 11, times the intensity of the collimated 

p 

1 ight beam used for illumination, IH may be written as 

11,2 = a 0 + a l I + a 2 j2 (13) 

where the coefficients a Q , a^, and a 3 include all the parameters which 
appear in Eq. 12 and the intensity of the light beam used for the illumination. 

B. The Mo ire Image 

The moire image is evaluated by multiplying the intensity transmittance 
distribution of the transparency and that of the grating, and averaging the 
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product over a pitch distance, namely 0 

5+1 

M(r,A,R,# ) = / m 2 G d£ (14) 

5 

o 

where ITI is given in Eq. 13 and G, the intensity transmittance distributi 
of the grid is represented by a Fourier series 7 


on 


G(r v ] 1 +1 y s in[ 2 n( 2 n+l ) (C+x/P+ye/P )] 
- 7 m L. £n+I 


* n=0 


(15) 


where x/P and ye/p represent the phase shifts related to the relative 
translation and rotation of the Ronchi grating. G is not affected by 
diffraction since the licjit intensity is measured directly behind the grating. 

When Eqs. (11), (13) and (15) are substituted in Eq. (14) and the 
integration is performed, we get: 


M(r , a,R,0 ) = A + B l — [ irU .. A *( 2n ^l ) J + 

n=0 (2n+l y 


4 c ? cos(tiua*[ (2n+l ) 2 - (2m + l ) 2 ]) cosUua! ( 2k+l ) 2 -(2£+l ) 2 ]) , , 

m # n,k,*-0 (2n+l)(2m+l) (2k+l)(2A+l) «(m,n,k, 5 t) 


( M 


“ cos[2*(2n+l)(£-+ + JL)] 

+ D I 2-5 — B E cos[nU 4 *( 2 n+l) 2 ] 

n =0 (2n+ir 


O 


(16) 


+ E l 

m,n,k ,£=0 


cos[2.(2n+l) (J + f - * 5*-)] cos[«A*(2ni*l) 2 ] 
f2n+T5 [2m+T5 




2 


a n a i u 3a ? u 2 a 9 u 2 „ 2a.u a 0 u 2 4a 0 u 

where A * ^ + -j jf-, B = C = D = -j- + E = , and 

6 (m,n ,k ,i ) = 6[ (n-m)-(k-£ )] + 6[ (n-m)-(A-k )] + d[ (n+m)-(£+k )] - <s[ (n-m)-(k+£+l )] 

- 6[ (m-n)-(k+£+l )] - fi[ (n+m)- (k-£-l )] - 6[ (n+m)- (£-k-l )] , where each 6 on the 
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right hand side of the equation is the Kronecker 6. 

As seen from Eq. 16, the moire image is strongly affected by diffraction 

(associated with the first Ronchi grating) and by nonlinear recording. Since 

the conventional deferred def lectometry suggests analyzing the moire pattern 

ill? 

by measuring the transmitted i ntensity x ±,± , it is interesting to display 
the moire image as calculated in Eq. 16. Figures 5a-51 illustrate the moire 
intensity vs the phase (x/p + ye/p), for various grating-photograpni c plate 
distances. All curves shown are without a phase object i.e. u = 1 and 
it> x = 0 . The distances are a = a + fA^ where f is a fraction that was 
varied in the range of 0.00-0.50. a^ and a^ are defined in Eq. 2. In each 
figure the two upper curves are the results for Agfa-Gevaert plates, prefogged 

O 

with E q = 12 erg/cm and exposed, through the grating, by a collimated beam 

2 2 

with E ’ =8 erg/cm and E' = 18 erg/cm for the solid and dashed curves, 

respectively. The lower solid curves are for Kodak plates with E Q = 0 and 

2 2 

E' = 55 uJ/cm . Exposure of E' =8 erg/cm falls in the linear region of 
the Agfa-Gevaert |T f - E curve, thus the E‘ = 8 erg/crr£ moire intensity 
curves exhibit only effects of diffraction. The effects of nonlinear 

O 

recording may be observed by comparing the E ' =8 erg/cm curves with the 
E' = 18 erg/cn£ and with the E‘ = 55 uJ/crn^ curves which are affected both 
by diffraction and nonlinear recording. 

For f = 0 a triangular wave form is obtained, both for the linear and 
nonlinear curves (Fig. 5a). At this distance, the exact dependence of the 
moire image on the phase is known, thus neither diffraction nor nonlinearity 
affect the accuracy and sensitivity of the system^. As f increases, the 
behaviour is changed: the maxima and minima are rounded, the contrast is 
reduced and higher spatial frequency waves are superimposed on the triangular 
wave. These changes are more pronounced in the nonlinear curves, and thus may 



- 10 - 


introduce hicjier uncertainties in the deflection measurements. At f = 0.25 
the wave form is again triangular, but with lower contrast as compared to that 
of f * 0. At f = 0.50 the intensity is independent of the phase and the 
contrast is zero. These results are of practical importance since 
experimentally it is hard to adjust the spacing between the Ronchi grating and 
the photographic plate to be exactly one of the Fourier image planes, 
especially in the case of an arbitrary phase object where such a single 
spacing, for the whole object, does not exist^’^. 


C. Heterodyne Readout 

As in the case of the on-line heterodyne moire, the heterodyne readout is 
accomplished by translating the Ronchi grating relative to the photographic 
transparency in the x-direction. The motion will cause the moire fringes to 
move in the y-direction. Assuming a linear motion of the grating with speed 
V, the phase of the moire image will vary with time as 


♦o-o - 2 ’ ( f * * nr + h 


(17) 


where Q = 2v — . If the h idler harmonics of the signal given by Eq . 16 
P 

are filtered out, the signal is given by: 


M T1 (r,R,A,t,tf x ) = cos[2r (jj + + (DcosUua*) + 


(18) 


+ f ? cosOuA*(2ninr] cos{hua*[ (2k + 1 ) -Mj ]} , k )} 

E ten 2m+l (2k+l ) (2t+l ) 

m,K ,^=u 


The fundamental frequency component of Eq . (16), for the reference 


transparency is: 
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Mpi(y,A*,t) = C0S[2n(~' + ^ + %Z t)](DC0S(irA*) + 

2 /On -Li ^2- 


R1 
+ E 


P P 2 it 

\2 


( 19 ) 


y c ?s[itA*( 2 m+l ) ] cos{ 1 A*L[ 2 k^l)^( 2 A , n) £ ]} , k £)} 

fe „_ n 2 m+l ( 2 k+l)( 2 £+l) 6 im,n-o,k ,£ )i 


m,k,£=0 


Comparing the transmitted signals with and without the phase object, it is 
seen that the effect of refraction is to introduce a phase shift trt(r) given 
by 


(Mr) a 
A¥(r) = 2 * 


( 20 ) 


already discussed in connection with Eq . 10. The procedure for measuring this 
phase shift is similar to that described for the on-line heterodyne 

q 10 

defiectometry- 5 . The phase of the signal, relative to a reference signal 

is measured twice at each point on the combination transparency-grating: once 
with the transparency recorded with the phase object and once with the 
transparency recorded without the phase object. The difference between the 
two readings gives the absolute phase shift given by Eq . 20. The reference 
signal is the output of a photodetector placed at a fixed point in part of the 
fringe pattern not affected by the phase object. 


D. Advantages of the Heterodyne Readout 

The theoretical analysis above shows some of the advantages of the 
deferred electronic heterodyne readout technique over the conventional 
deferred intensity readout method: 

1. The errors introduced because of inaccurate positioning of the 

photographic plate precisely on one of the Fourier planes (diffraction and 
the nonlinear recording effects) are small compared with those in 
intensity measurements. This is true since the electronic phase recording 
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of the fundamental frequency signal is insensitive to small fringe 
contrast variations and to the fine structure of the fringes (see Fig. 5) » 

2. Variations of intensity over the fringe pattern due to optical 
imperfections and other effects do not affect the phase measurements since 
the phase is computed as a quotient of two intensity terms 9 . However, 
for intensity measurements, this effect may lead to significant errors. 

3. The electronic heterodyne method has the ability to interpolate 
continuously between fringe maxima, and has computer compatibility. 

4. In principle, by using the third harmonic term of Eq. 16 (with the time 
varying phase term 3 ^ t) the electronic phase of 3ai l> will be measured, 
which is equivalent to tripling the sensitivity of the system 9 . 


E. Spatial ana Angular Resolution 

Diffraction effects, detector aperture and electronic instrumentation 
determine the spatial and angular resolution of the system. Since the 
nonlinear photographic recordings do not affect the resolution of the system, 
the results of the spatial and angular resolution obtained for the on-line 
def lectometry are valid for the present deferred case 9 . 

The spatial resolution of the system is 

^ = p* + 2a* (21) 


where p* = p/p is the dimensionless size of the detector. 
The angular resolution is given by: 

U6) min = 3WT* 


( 22 ) 


where x* = x/p. 
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F. Some Practical Consioerations 

It is important to point out that the final accuracy and sensitivity of 
measuring phase ^ is strongly dependent on the manner in which it is 
practically performed. 

1. The electronic processing of the photographic transparencies is performed 
on a separate bench. The procedure for measuring /up (Eq. 20) is to record 
the electronic phase, relative to that of a reference signal twice: once 
when the phase object transparency is placed in front of the grating and 
once with the reference transparency in front of the grating, is then 
obtained by subtraction of the readings. It is clear that the accuracy of 
the system and the ability to measure small deflections is dependent on 
the capability to accurately place the two transparencies precisely in the 
same position relative to the grating and to the detectors. Small 
misalignments will result in considerable errors in the measured phases 
and thus in fringe deflections. 

2. Practically the phase meter output is dependent on the visibility of the 
fringes and on the distortions of the fringe profile caused by the 
nonlinear recording. Increasing the visibility, on the one hand, reduces 
the error of phase reading. On the other, hand it increases the effects of 
nonlinear recording which affects the measurement as well. Thus in 
forming the transparencies, one must compromise between those two effects. 
In the previous section it was mentioned that by using the third harmonic 

component of the moire signal, it is possible to tripple the sensitivity of 
the system. It is worthwhile to point out that practically it is hard to 
materialize this useful feature since the third harmonic signal is weak and 
thus the s i gnal-to-noise ratio is low. 
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III. EXPERIMENTA L 

The deferred electronic heterodyne technique was evaluated by measuring 
small fringe shifts caused by a weak phase object. The phase object used was 
a large crystal of KD*P whose refractive index distribution could be changed 
when voltages are applied to the electrodes mounted on the crystal. A 
detailed description of this phase object, a so-called gas-flow simulator, is 
given in Ref. 15. The flow simulator was chosen to be tested because of two 
reasons: a. The flow simulator was previously investigated by the electronic 
heterodyne holographic interferometry and by the on-line electronic heterodyne 
moire def lectometry methods^®. The availability of these measurements, 
which serve as a comparison makes it possible to evaluate the system, b. The 
flow simulator is a very weak phase object which exhibits a fine structure 
close to the limit of the angular and spatial resolution of the system. The 
ability to measure such a phase object validates the theoretical predicted 
performance . 

The moire setup for recording the phase object was the same as that 

Q 

described previously for the on-line experiment except that the second 
grating was replaced by a photographic plate. We used Agfa-Gevaert 1QE75 
plates. The phase object transparency was recorded when the high voltage of 5 
kv was applied to the crystal. The reference recording was performed with the 
crystal in place, but without the voltage. 

Other parameters in the experiment were the same as those in Ref. y, 
namely: the separation between the grating and the photographic plate was a = 
32 cm, the grating pitch was p = 0.17 mm thus x* = 3.27x10“^, a* = 7, the 
spatial resolution, calculated from Eq. 21 was ax = 3.1 mm and the angular 
resolution, from Eq. 22 was a n = 1.5xl0 -b rad. 

The electronic processing of the developed photographic plates was 
performed on a separate Dench. The distorted moire fringes were observed when 
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the phase object transparency was placed in front of the Ronchi grating and 
illuminated by the collimated beam. The separation between the fringes, in 
the region away from the phase object measured p' = 4.47 mm. The measurement 
was as follows: while the grating was in motion and the reference detector was 
fixed at one position, the test detector was moved in the y-direction until 
the phase meter output covered exactly 360° indicating that the detector had 
moved a distance equal to p', the period of the moire fringes. 

The fringe motion was achieved by engaging the grating, as described in 
Ref. 10, to an electrical vibrator. The detectors and the rest of the 
electronic system were the same as those described in Ref. 10. The phases 
relative to the reference detector, were recorded by a the test detector along 
the center line of the crystal's image, at stations separated by 1 mm. 

At this point the phase object transparency was replaced by the reference 
transparency and phases were recorded exactly at the same stations as before. 
The exact positions of the detector were determined by an interferometer. A 
point by point subtraction of the phases, with and without the phase object, 
yield the phase shift as a function of position and thus the index of 
refraction gradient vs. position. The results are shown in Fig. 6a. The 
results show that the measurements were performed close to the angular and 
spatial resolution limits. Indices of refraction of the order 2xl0 - ^ cm“^ 
were measured and the peak structure of spatial scale of 3-7 mm was resolved. 

Figure 6b displays the results measured with the on-line heterodyne 
def lectometry method. The results obtained with the heterodyne holographic 
interferometry method are shown for comparison in Fig. 7. The figures show 
that within the error limits, the results obtained by the three methods agree 
very well. It is seen that the deferred results are consistently higher, 
except at the top of the third peak, than those obtained by the two other 
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methods. The source of this error may result from the transparency 
misalignment which could have caused a constant phase shift over the whole 
field. 

Reproducibility of the measurements was then tested on a few photographic 
plates taken of the same object with different exposure levels. We found that 
the phase was only reproducible to within -*4* over the field. These errors 
may be attributed to mechanical misalignment of the plates and to contrast 
variations and non-linear recording effects. More experimental work has to be 
performed in order to investigate these effects. 

For sinusoidal grating motion, rather than motion at a constant velocity, 
the output of the electronic phase meter does not vary linearly with 
refraction angle, thus calibration of the system is required 9 . In the 
present experiment the calibration procedure was the same as that followed in 
the on-line experiment, except that it was performed with the reference 
transparency in place (Fig. 3) rather than with the actual KD*P crystral . 
Similar calibration curves to those shown in Fig. 7 of Ref. 9 were obtained. 

Since the accuracy and sensitivity of the system is largely dependent on 
the transparencies alignment, a special effort was made to put the reference 
transparency precisely in the same position, relative to the Ronchi grating 
and to the detection system, as that of the phase object transparency. This 
was achieved by the following: a. The transparencies were mounted on a plate 
holder, fixed to the optical bench. The use of the plate holder ensures that 
the surface coordinates of both transparencies was the same. b. The rotation 
angle between the reference transparency and the grating was adjusted by 
measuring p', the separation distance between the moire fringes. The angle 
was varied until the value of p' reached the same value as that measured with 
the phase object transparency, namely p'=4.47 mm. c. A geometrical point, 
photographed on both transparencies, was chosen as a reference point for the 
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alignment procedure. Before removing the phase object transparency, the x, y 
coordinates of this point were measured. In order to align the transversal 
position of the reference grating, it was transversal ly moved until the x, y 
coordinates of the reference point attained equal values to those measured on 
the phase object transparency. 


IV. CONCLUSIONS 

It has been shown that the electronic heterodyne technique can be applied 
for accurate and sensitive readout of moire fringes of unsteady phase 
objects. Theoretically it was shown that the accuracy and sensitivity of the 
system are weakly affected by diffraction and non-linear photographic emulsion 
characteristics. The analysis showed that non-linear recording does not 
affect the spatial resolution of the system. The spatial resolution is mainly 
affected by diffraction. 

The system was evaluated experimentally by mapping a weak phase object. 

The results compared favorably with observations performed with the electronic 
heterodyne holographic interferometry and with the on-line heterodyne moire. 
The theoretical analysis and the experimental measurements have shown that the 
performance of the the heterodyne deferred moire is of the same order as the 
on-line heterodyne moire. 

The method offers many advantages: low stability requirements, adjustable 
sensitivity, not affected by noise sources, low cost and conputer 
compatibility. To conclude it can be said that the method is suitable and 
attractive to use for non-steady phase objects. 

The deferred electronic heterodyne moire requires continuing research and 
development. More experimental investigation is necessary in order to better 
understand the reasons for the low reproducibility observed. It is suggested 
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to consistently study the effects of photographic plates exposure levels and 
developing processes on the system performance. In the present research the 
only photographic emulsion used was the Agfa-Gevaert 1QE75 one, not 
necessarily the most suitable for moire. Other emulsions with different 
characteristics should be tested as well . The errors caused by rrechanical 
misalignment of the reference transparency when placed in the position of the 
phase object transparency, have to be reduced. This can be achieved by 
including better alignment facilities on the processing bench. Also, as in 
the on-line case, the sinusoidal motion of the Ronchi grating is undesirable 
because of the need for calibration. A technique for linear motion has to be 
developed. 
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PART 2: MEASUREMENTS OF 3-D DENSITY FIELDS WITH OPAQUE OBJECTS 

BY MOIRE DEFLECTOMETRY 

I. INTRODUCTION 

Tomographic reconstruction of 3-D refractive-index fields from measured 
ray deflections by moire method was demonstrated in Ref. 16. It was shown 
that if the data on deflections is complete i.e. the field is not blocked by 
an opaque object and deflections of 360° coverage are available, the inverse 

17 10 

Radon transform * may readily be used for field reconstruction. The 
advantage of the def lectometry method is that the measured deflections are 
directly used in the inversion transformation. This is not the case with 
interferometry. Interf erometri c measurements, namely Jndx where n is the 
index of refraction and x is the distance coordinate along the ray, should be 
mathematically processed in order to obtain numerical derivatives. This 
process may add errors to the reconstructed field. 

In some cases it is desired to use moire def lectometry for measuring 3-D 
density fields where an opaque object is present in the field. This situation 
occurs when def lectometry is used to study density fields in compressible gas 
flowing around a test object or to study the temperature fields around heated 
opaque objects, etc. In such cases the object blocks many of the optical rays 
used for the measurements and thus part of the data is missing. Since the 
missing data is necessary for evaluating the inverse Radon transform, the 
direct use of the reconstruction technique is impossible. 

The same problem arises when interferometric incomplete data is used for 

1 9 20 

field reconstruction; Zien et al. , Vest and Prikryl and Medoff et 
? 1 

al. proposed to solve the problem by using algorithms based on iterations. 
For reconstruction it is necessary to extrapolate the projection data into the 
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region of missing data. The artificial field evolves during the iterative 
process by transforming repeatedly betveen the image domain and the projection 
domain, applying measured information and other criterions at each iteration. 
Empirical studies of this algorithm shows 20 that the algorithm does not 
converge. The studies indicate that the field reaches a best estimate after a 
number of iterations and then diverges. A similar iteration technique was 
previously used with deflections data to reconstruct the flow field around a 
cone at an angle of attack . The problem of algorithm divergence was also 
encountered in that study. 

In many practical situations the opaque object presented in the field is 

of circular or elliptical cross-section. This occurs in aerodynamics of 

circular bodies, like cones, at an angle of attack, or in heat-transfer of 

circular pipes, etc. For such bodies it is shown, that by using the inverse 
18 23 

Cormack transform ’ the field can be reconstructed directly without the 
need for iterations. The Cormack transform is based on circular harmonic 
decomposition of both the image space and the projection space and solving the 
relation between these deconpositions. It is interesting to note that Cormack 
transform is an extension of Abel transform (for circular symmetric fields 
Cormack transform reduces to Abel transform). 

A brief description of the inverse Cormack transform is given in Sec. II. 
It is shown how this transform is used with deflections data. An algorithm, 
based on Cormack transform, is developed and evaluated in Sec. III. The 
Appendix contains listing of the computer program. 


I I. THEORY 

Cormack's approach is based on the fact that for each r the object 
function f(r,e) is periodic in e with period 2tt, where r and e are the polar 
coordinates defined in Fig. 8. Therefore, f (r,e) can be expanded in Fourier 
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Series of the form of 


f (r,e) = l f„(r)e 


i He 


( 1 ) 


l = -°° 


where 


f £ (r) = 27 / f(r,o)e~ 1£6 de . 


( 2 ) 


Similarly, the interferometric projection x(p,b) can be expanded as 


x(p,b) = l x (p)e 

£= -<=° 


ub 


(3) 


with 


, rZ ir -Ud 

(p) = T - x (p , jb) e db 


Y p ' = ^ o 


(4) 


x(p,b) is defined by the line integral 


00 2 IT 

x(p,b) = / / f (r,e) a[p-rcos(&-b)]rdrde 
0 0 


(5) 


where the integration is performed over the infinite plane but the delta 
function in this equation has the effect of reducing the area integral to a 
line integral along p=rcos(&-d) , shown in Fig. 8. <t> is the direction of 

integration and p is the distance of the line of integration from the origin. 

Cormack transform gives the relationship between the Fourier coefficients 
x^(p) and f £ (r) 18 : 


x £ (p) 


~ V r > Vr> 

2 / i - „— I- rdr 


( 6 ) 
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where T (=•) is the Tschebycheff polynomial of the first kind of order z defined 
by: 


T*(x) = 1 


COS(£COS *x) 
-1 


|x| < 1 


coshUcosh x) |x| > 1 


(7) 


Note that if f(r,e) is rotational symmetric, f^(r) is non-zero only for 
£=0 and in this case Eq. 6 reduces to the Abel transform, since T 0 (x)=l for 
al 1 x . 

For the reconstruction of the object from its projections, we must 
determine the inverse Gormack transform i.e. to express f ^ ( r ) explicitly. The 


solution, in one of its forms, is given by^: 


dX £ (p) D 

f t (r) - 

r \I7^ 


dp 


( 8 ) 


Interesting features of this result are: 

1. Only values of x(p,<6) for p>r are required to reconstruct f(r,©), thus Eq. 8 

can directly be used for reconstructing fields around circular objects 

since the data gap in the region p<r is not required for the integration, 
dx (p) 


2. Tne derivatives 


dp 


-, which appear in Eq. 8, are exactly the g, ' th 


Fourier coefficients of the function = x'(p,d). The proof is straight- 

9p 

forward, since x(p,d) is periodic in d, x'(p,d) is also periodic, 
therefore can be expanded in Fourier series of the form of: 


x (p,6) = l x (p)e 


ut 


(9) 


£= — » 


where 


2 TT 

x (p) = / x (p,/t>)e ^ d d<!) 

0 


( 10 ) 
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by applying Eq. 4, Eq. 10 may be written as 

. , ^ 2 ir . . dx (p) 

\ (P)= ^3P 0 ^ X( P’^ )e 6t) = '~W~ (11) 

It turns out that 3X fP v -- is measured directly by def lectanetry 1 and 

dP 

therefore no differentiation of experimental data is required prior to the 
expansion of the measured data in Fourier series. 


III. COMPUTER SIMULATION 

The reconstruction procedure discussed in the preceding section was 
evaluated by a computer simulation. The phase object studied was analytically 
defined and the algorithm was used to reconstruct the object from conputed- 
generated deflections data. The object function was zero everywhere outside a 
circle of radius R and an opaque object of circular cross section of 
radius R was placed in the object field. In order to have a basis for 
conparison of the accuracy of the reconstruction obtained in the present work 
with that obtained by using an iteractive method, we chose to study the same 
object as that studied by Vest in Ref. 20. 

The phase object studied, shown in Fig. 9 is given by: 


f(x,y) = exp ( = fl (x - i^&& ) +0.5 exp C 6 K** 0 ; 6 I V 3 } 
l-(x 2 +y 2 ) l-(x 2+ y 2 ) 

+ exp^Lx^Wl) ♦ o.5 exp f^Vo-e) 2 ] ) (i 2 ) 

l-(x 2+ y 2 ) l-(x 2+ y 2 ) 


for 


R < r < R 
c - o 


f(x,y) = 0 for 0 < r < R and r > R . 


and 
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The conputer simulation program, listed in the Appendix, consists of the 
following steps : 

1. Numerical generation of deflections data, x'(p,£). For this purpose 
f(x,y) was expressed as a function of p, s, and <t> according to the 
coordinate transformation 


x = pcosd - ssinb 
y = psi nd + scosi 


(13) 


where s is the coordinate along the ray direction - see Fig. 8. Next, the 
object function was differentiated with respect to p and numerically 
integrated along s, yielding x ' (p • jb) . 

It was assumed that we have N projections spaced by equal intervals 

A 

U = 411 . jhe number of data points along the p axis is M thus the 

R -R 

increment in p is ap = — jvjZ - ! - • 

2. For each p, x ' (p , 6) is expanded in Fourier Series in the azimuthal 
direction (Eqs. 3 and 4). The number of terms retained is £ m ^=L. 

3. Equation 8 is used to calculate f^r), the Fourier coefficients of the 
azimuthal expansion of the object function. In Eq. 8 the integrand 
becomes infinite as p->r. The singularity at p=r can be evaluated by 
dividing the integral in Eq. 8 into two parts^ 


dx^(p) 

I "oljlVfl 

TT j 


( 


dp = - 


2 2 
P -r 


dxJp) 

1 f r * h 

* r 


J2 2 ' 

1/ 


dp - 


I/O 

n r+h 


dx »( p ) n 

T (£) 

dp Vr' 

( 14 ) 


dp 


For small values of h, the first integral is approximated by 

l dx 0 ( r ) _i h 

- ^ [ — |p— T (1) cosh i (l + |)], thus Eq. 14 becomes 
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dx (p) 

i R — h — 

1 , o dp 

" r IT 2 


Ld! i 


{7Z 


dp = 


1 d \(r) 

_ is — £ j 

it 1 op 


dx (p) 

1 u R ~ — T (-B-) 
(l)cosh -1 (l+ $) * r 0 -SB i. L 




/7^ 


"i 


(15) 


dp} 


The approximation error is of the order of O(h^) . 

4. The object field function is calculated from Eq. 1 using the computed 

V r) - 

A summary of the reconstruction results is given in Table 1. 


TABLE 1 Reconstruction Errors Obtained by Inverse Cormack Transform 
and Iterative Convolution Methods* 


Opaque object 

N 

M 

L 

Number of 
grid points 

Maximum 
Error, % 

Averaged 
Error, % 


35 

29 

8 

1015 

182.0 

9.5 


30 

25 

6 

750 

15.8 

1.8 

Circle ( R =0.3R ) 

30 

25 

5 

750 

8.7 

1.6 

C 0 

30 

30 

5 

900 

8.5 

1.7 


Iterative 

convolution 

19.5 

1.8 


35 

11 

10 

385 

22.0 

2.7 


35 

16 

6 

560 

6.5 

1.3 


35 

16 

8 

560 

5.4 

1.3 

Circle (R =0.6R ) 

35 

16 

10 

560 

2.0 

0.3 

C 0 

30 

16 

10 

480 

5.3 

0.6 


30 

16 

5 

480 

8.5 

1.5 


Iterative 

convo lution 

27.5 

3.7 


iterative convolution results are from Ref. 20. 


In the table the maximum error is defined as the maximum value of the 
differences between computed and actual values of the object function at each 
point, divided by the maximum value of the function that occurs in the object 
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domain. The average error is defined as the average of the absolute values of 
the errors over the whole grid. For comparison, the errors calculated by the 
iterative convolution method are shown in the Table as well. From the Table, 
it can be seen that for this example, the inverse Cormack transform is 
considerably more accurate than the iteration method; maximum errors for 
R c =0.3R q and for R c =0.6R o are 8.5% and 2. 0 % respectively as compared 
to errors of 19.5% and 27 .5% obtained from the iterative convolutin method, 
respectively. 

I V. CONCLUSIONS 

The series expansion approach used in our investigation yielded a fairly 
efficient reconstruction code. It was shown that the use of the inverse 
Cormack transform for reconstruction of fields with circular missing 
deflections data is quite accurate and convenient for use. The accuracy of 
the algorithm was tested by applying it to data generated by numerical 
simulation. The relatively high accuracy found is mainly due to the fact that 
no iterations are necessary. 

Our study was limited in this project to asymmetric fields around circular 
cross sectional opaque objects. It is important to continue and expand the 
algorithm for elliptical cross-section opaque bodies (which often occur in 
aeronautics, when circular bodies are at an angle of attack) and for general 
shaped bodies. Also more analytical studies have to be done in relation to 
the efficiency of the algorithm. 

It has to be emphasized that in the present study we have not made any 
effort to reduce conputational time. The principal factor in determining 
computational time is the number of grid points and the number of terms 
retained in the Fourier series. Thus by using the Fast Fourier Transform 
technique the computational time may be considerably reduced. 
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APPEfOI X 

LISTING OF THE COMPUTER RECONSTRUCTION PROGRAM 


KF = Number of Fourier coefficients in the series expansion. 

NP = Number of data points along p axis. 

NFI = Number of projections. 

RCRO = R c / R 0 • 

DFI = The interval in projection angle. 

DP = Increment in p. 


MAIN PROGRAM 

THIS program reconstructs density fields around circular 

OPAQE OBJECTS FROM DEFLECTIONS DATA. THIS PROGRAM EMPLOYS 
THE INVERSE CORMACK TRANSFORM 


IMPLICIT REALS81 A-H,0-Z) 

INTEGER IER 

EXTERNAL F , RFF , A I FF , FOU , FOUR 

DIMENSION P ( 30 ) , FI ( 60 ) ,RDN ( 6 0 > 30 ) , RLMD ( 30 ) , A I L MD ( 3C ) , RD ( 60 ) 
DIMENSION FID( 60 ) ,RAF ( 30 , 60) , RI AF ( 30 , 60 ) , RREAL (60,30) ,RAIM(60,3 
DIMENSION X X ( 60 ) ,YY(60) ,FX ( 4 ) ,FY< 4) ,C( 4, 3 l 
DIMENSION RF ( 60 , 30 ) , A IF ( 60 , 30 ) 

COMMON/FM/SP , SN , CS , PM , PCS, PSN 
COMMON/RFM1 /RLMD , A I L MD , NP 
COMMQN/RFM2/P , R , RSQ , I R 
COMM ON/ A A A/RD ,FI , NF I P 
COMMON/ B3B/LR , NF I 

DEFINITION OF THE DENSITY FUNCTIONTO BE RECONSTRUCTED 


A . 4 . . 4 . .4. . 4 - . 4 . . 4 . ... 

WV '<* V V V V V -S” W V V V V 


.4. ,4. .1, .1, .4. 

•i” t r vv v 


■. V. .A. *.4 - *4. wAw .4. .4. .A. 

V W V V Y“ V V V T W V V V V V 


PE XP ( X, Y) =DEXP{ -6.« ( ( X-0 . 6 ) Z*2 + Y**2 ) / ( 1 * - ( X*X+ Y* Y ) ) ) +0 . 5SDEXP< - 
0 ( ( X*0. 6) «*2+Y*S2 ) / { i . - C X*X*Y*Y ) ) ) +0EX°(-6.4( X=C=S2* ( Y-0 . 6) X=*2) / ( i 
#XS^2 + Y=S*2 ) ) ) +0. 5X=DEXP (-6.* < ( Y + 0.6 ) /( l . - ( X* X+YSY ) ) ) 


CALL UERSET ( 0 , LEVOLD ) 
PYE=3. 141 50265358979 
CR=6. 0/1 74.0 
CH=0.01 
AERR = 0 . 

F AERR=0 • 

RO= 1.00 
KF= 1 0 
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RCR0=.6 

RERR=.01 

FRERR=0.01 

NP=I 1 

SRO=RO*#2 

RC=RCRO*RO 

NFI =25 

NFIP=NFI+1 

TTL=NFI*NP 

OFI=2.*PYE/'NFl 

OF I D = 360. /NFI 

NPH- NP— 1 

DP=(RO-RC)/NPM 

DO 1 0 M= 1 , NP 

AHM=M— 1 

P(M)=RC+DP*AMM 
10 CONTINUE 

DO 11 N=1,NFIP 
AMM=N- 1 
FI|N)=DF1*AMM 

F I D( N } =OF I OvAMM 

11 CONTINUE 

WRI TE { 5 | 70 J NFI ,NP,KF ,CH, TTL 

70 FORMAT (2X, * NF I , , NP , , KF , , CH , , TTL * , 3 { 3X , I 2 ) , 3X , 2F 1 O . 2 ) 


CALCULATION OF THE DEFLECTION PROJECTIONS (INTEGRAL) 
OF THE GIVEN FUNCTION PEXP(X,Y) FOR THE DEFINED 
VALUES FI AND P WHICH WE CALL RDN(PHI,P). 


13 

12 


19 


53 

50 
54 

51 


WRI TE ( 6 | 5 1 ) 

DO 12 N=1 , NF I 
CS=DCOS(FI IN) ) 
SN=DSIN(FI ( N) ) 
DO 13 M=1 , NPM 
PM=P(M) 

SP=P( «|t$2 
PCS=P( M)SCS 
PSN = P ( M ) 4=SN 
B=DSQRT( SRO-SP) 
RA=-8 


RDN(N,M)=DCADRE(F,RA,3, AERR ,RERR, ERROR , IER) 

CONTINUE 

RDN ( N | NP ) = O « 0 

CONTINUE 

DO 1 9 M= 1 , NP 

RDN ( NF IPjM) =RDN ( 1 , M) 

CONTINUE 


WR I TE ( 6 » 50 ) (FID(N) ,N= 1 ,NFIP) 
WRI TE ( 6,50 ) ( P( M) ,M= 1 ,NP) 


DO 53 M=1 , NP 

WRITE (6,54) ( RDN( N,M) ,N=l , NFI P ) 

CONTINUE 


FORMAT ( 3X, l OF 10. 4// ) 
F ORMAT ( 3X , 1 OF 1 0 . 4 ) 
FORMAT! 3 X , *PHI = 


P = 


• ) 


RDN = 
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expansion OF RDN(FI,P) IN FOURIER SERIES 
WITH P ASA PARAMETER 

$ 4= * $ ❖ ❖ * * * * * S 4 = 4= ^ # 


» ^ .•* A ^ « 


' * 1 * W T V '(* "V* 'i* V 


250 


26 
I 17 


57 


55 

56 
68 


I 

I 


i 


220 


210 

200 


SS-0 - 

TT=2 ,*PYE 

do 17 m=i,npm •; 

DO 250 N=1,NFIP OF IROO^ -H l'\l JTY 

RD{N)=RDN( N,M) 

CONTINUE 
OO 26 LR=1 ,KF 

RRE AL ( LR , M ) = DCAORE ( FOU , SS , TT , F AERR * F RERR, ERROR , I E R ) 
RREAL( LR * M ) =RREAL ( LR , M ) /PYE 

RAI M(LR,M)=DCADRE( FOUR, SS, TT,FAERR,FRERR, ERROR, IER) 

RAIM( LR,M) =RAIM(LR,M ) /PYE 

CONTINUE 

CONTINUE 

OO 57 LR=l,KF 

RAIM(LR,NP)=0.0 

RREAL ( LR , NP ) = 0» 0 

CONTINUE 

00 55 M- 1 ,NP 

WR I TE ( 6 , 56 ) M,P(M) 

WRITE (6,54) ( RREAL (LR,M) ,LR=1,KF) 

WRITE (6, 54) (RAIM(LR,M) ,LR=1 ,KF) 

CONTINUE 

FORMAT(4X, 'FOURIER COEF. OF RDON , M= • , I 3 , 3X , F6 . 3 ) 
FORMAT ( 3X,F14.6,2X,F14.6,2X,I4,2X,I4) 


CALCULATION OF THE FOURIER SERIES EXPANSION OF THE 
DENSITY FUNCTION BY USING THE INVERSE CORMACK TRANSFORM 

s> 4: « s£ S * * # « 4 ***$<:$**:$* 4 4 4 4 « 4 « S « S <5 « 


H=CH«RC 

DO 200 I R= I , NPM 
R=P( IR) 

RSQ=R#«2 
AR=R+H 
ARH= AR/R 

ACSH=DLOG( ARH + DSQRT ( ARHSS2- 1 ) ) 

DO 210 LR= 1 , KF 
OO 220 M=I,NP 
RLMD ( M ) =RRE AL ( LR , M ) 

AILMD( M) =RAIM(LR,M) 

CONTINUE 

RF ( LR , I R ) =DCADRE l RFF , AR , RO , AERR , RERR , ERROR , I ER ) 
AIF< LR , IR) =DCADRE ( A IFF, AR, RO , AERR , RERR , ERROR , IER) 
RF ( LR , IR)=(-1/PYE)«(RLMD.( I R )$AC3H*RF (LR , IR) ) 
AIF(LR , IR) = (- 1/PYE) *{ AILMD ( IR ) SACSH4-A IF ( LR , IR ) ) 
CONTINUE 
CONTINUE 


RECONSTITUTION OFTHE DENSITY FUNCTION FROM THE 
CALCULATED FOURIER COEFFICIENTS 


«A» 2 ** "4 »•* «4 A »•» «>» «*• A *•» •*« ■'» <4 «4 ••• A -A> <4 «A» *4 i 4 -J 4 MU MW 

w a* v v-y v ir v *r v V V v*r v V v* -c v v w '«* w w*c w w w* *<» -v* v y y 




APMR=0.0 

DO 400 I R= 1 , NPM 

DO 312 N= 1 , NF I 

RAF ( 1 R,N ) =C.5*RF ( 1 , IR ) 

RIAF( IR,N)=0. 

OO 300 L=2,KF 
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OF PX 


IS 


QUALITY 


300 

444 

445 

18 

312 

400 

987 


Z=L-1 

RAF{ I Rj N) = R AF ( I R , N ) + RF ( L , I R ) *DC0S < Z4F I (N) )+AIF(L, I R ) *DS I N ( Z4FI (N) ) 
CONTINUE 

X=P( IR) ^OCOSC FI (N) ) 

Y=P< I R ) 4=DS I N I F I (N) ) 

IF(DABS( l.-X4*2-Y4:$2) .LE.CR) GOTO 444 
PP=PE XP ( X » Y ) 

GO TO 445 

PP=0.0 

CONTINUE 

PPMR=(PP-RAF(IR,N) )4=100. 

APMR=APMR+DABS ( PPMR ) 

WRI TE(6|18) PP,RAF( IR,N) ,PPMR,IR,N 

F 0RMAT(3X,F14.6,2X,F14.6,2X,F14.6,2X, I4,2X, 14) 

CONTINUE 

CONTINUE 

APMR= APMR/TTL , 

WRITE (6,987) APMR , TTL 

FORMAT ( 5X///’ AVERAGE ERROR= ' , F 1 0 . 4 , » NUMBER OF POINTS =',F5.1) 

STOP 

END 




wA* -JC Wi* 


4* 


FUNCTION F ( S ) 

IMPLICIT REAL*8(A-H,0-Z) 

DIMENSION P(30) ,RLMD( 30 ) ,AILMD(30) 
COMMON/RFM 1 /RLMD , A I LMD , NP 
COMMON/RFM2/P , R , RSQ , IR 
COMMON/FM/SP , SN , CS , PM , PCS , PSN 
COMMON/BBB/LR , NF I 
T=SP+S4$2 
Cj= 1 .-T 

IF( Q.LE. 1 .00-75 ) GO TO 111 

SQ=Q*42 

SCS=S*CS 

SSN=S*SN 

EF 1 = ( T+0.36-1 • 2=C ; ( PCS— SSN ) ) /Q 
IF(DABS(EF1 ) .GT.29.0) SO TO 10 
F1=DEXP(— 6.0SEF1 ) 

DF1 = — 6. 0*F1*( ( 2.SPM-1 .2SCS ) +2.$PM*EF1 ) /Q 
GO TO 11 

10 F 1=0 • 0 

DF 1 = 0.0 

11 EF2=(T+0.3b+1.2*( PCS- SSN ) ) /Q 
IF( DA8S( EF2) .GT.29.0 ) GO TO 20 
F2=0.5^DEXP (-6. 0$EF2) 

DF2 = - 6 . 0iF2i ( ( 2 .$PM+ 1 . 2*CS ) * 2 . *PM=5=EF2 ) /Q 
GO TO 21 


20 F2= 0 . 0 

DF2 = 0.0 

21 EF3=( T+0.36-1 .2*(PSN+SCS) )/Q 

IF (DABS(EF3) .GT.29.0) GO TO 30 
F3=DEXP(-6.0*EF3) 

DF3=-6 . 0£F3« ( (2.*PM-1 . 2*SN ) +2 ,4PM*EF3 ) /Q 
GO TO 31 


I 


- 31 - 


30 

31 

■ 

40 

41 

' 111 

222 

i 

i 

I 


F3= 0 • 0 
DF3=0.0 

EF4= (T+0.36+1 -24(PSN+SCS ) ) /Q 
IF ( DABS ( EF4 )»GT«29*0 ) GO TO 40 
F4= 0 . 54DEXP { -6 . 04EF4 ) 

DF4=-6.04F44{ (2.4PM+1 .24SNJ «-2.4PM*EF4 ) /Q 

GO TO 41 

F4=0.0 

OF4=0.O 

CONTINUE 

F=DF1 ♦OF2+OF3+OF4 .. . 

GO TO 222 
F=0 • 0 
CONTINUE 
RETURN 
END 

44 444*44444*44 $#44 444444 4444444 
FUNCTION FOU(X) 

IMPLICIT REAL48(A-H,0-Z» 

DIMENSION FI ( 60 ) ,RD(60) 

COMMON/ - AAA/RD, FI ,NFIP 

COMMON/BBB/LTR^NFI 

CALL FINTERCNFIP,FI ,RD,X,YYJ. 

LLL=LR-1 

FOU= YY4DCOS ( LLL4X ) 

RETURN 

END 

4444444444444444444444444444444 
FUNCTION FOUR { X } 

IMPLICIT REAL* 8 (A— H»Q— Z) 

DIMENSION FI160) ,RD(60) 

COMMON/AAA/RD, FI ,NFIP 

COMMON/BBB/LR > NF I 

CALL F INTER ( NFIP.jF IjRDjXjYY) 

LLL=LR-1 

FOUR=YY4DS I N l LLL4X ) 

RETURN 

ENO 

4444444444444444444444444444444 
FUNCTION RFF(X) 

IMPLICIT REAL48I A-H,0-Z) 

DIMENSION P ( 30 ) ,RLMD(30) ,AILMO(30) 

COMMON/RFM1 /RLMD, AILMD,NP 

COMMON/RFM2/P, R»RSQ, IR 

COMMON/BBO/LR > NF I 

LL=LR— 1 

MM=NP 

XR=X/R 

CALL FINTER(MM,P,PLMD,X,YY) 

RFF=YY4DC0SH(LL4DL0G(XR+DSQRT( XR442-1 ) ) ) /DSQRT ( X44 

RETURN 

END 



2-RSQ) 



FUNCTION AIFF(X) 

IMPLICIT REAL *8 { A-HjO-Z) 

DIMENSION P ( 30 ) ,RLMD{30) ,AILMD(30) 

COMMON/RFMl/RLMD, AILMD, NP 

COM MON/RFM 2/P , R , RSQ | IR 

COMMON/8BB/LR , NF I 

LL=LR— 1 

MM=NP 

XR=X/R 

CALL FINTERIMM, P, AILMD, X,YY) 

AI FF— YY*DCOSH{ LL^DLOG ( XR+DSQRT ( XR4:X=2— 1 ) ) ) /DSQRT ( X£S 2— RSQ ) 

RETURN 

END 


vWY-V'^'vvW'i' v v -v -y v 1 1* w -v “V v ' 


SUBROUTINE F I NTER ( N , XX , YY , X, Y ) 

IMPLICIT REALMS ( A-HjO-Z ) 

DIMENSION P(30) ,RLMD( 30) ,AILMD(30 ) 

DIMENSION FI(60I,RDN(60,30I 
COMMON/ A A A/RD N « F I 
COMMON/RFMl/RLMD, AILMD, NP 
COMMON/RFM2/P, R , RSQ, I R 

DIMENSION XXI N) , YY ( N ) ,FX (4) ,FY<4) ,C|4,3) 

NF=N 

IF! DABS I XX ( 1 )-X ) .LE.DABS ( 1 .D-6#XX t 1 ) ) ) X=XX{ I ) 

I F ( DABS I X— XX ( NF ) ) .LE.DABS I 1 . D-6*XX I NF ) ) ) X=XX(NF) 

IFIX.LT.XXI1 ) .OR.X.GT.XX(NF) ) GO TO 4 

OO 7 I=2,NF 

IF( X. GE. XX ( I ) ) GOTO 7 

K — 1 — 3 


M = 2 


IFII.EQ.2) K=0 
IFII.EQ.NF) K=NF-4 
I F ( I . E Q . 2 ) M=1 

IF ( I . EQ.NF ) M =3 
00 8 J= 1,4 
FX( J)=XX(K*J) 

FY( J)=YY(K+J) 

CALL IQHSCU(FX,FY,4,C,4, IER) 

Y=( <C( M, 3 )*{ X-FX ( M) ) +CI M, 2 ) I*|X-FX( M) ) +C ( M, 1 ) ) ❖ ( X— F X ( M ) ) +FY(M) 

GOTO 3 

CONTINUE 

Y=YY ( NF ) ; 

GO TO 3 

WRITE(6,5) NF, XX ( 1 ) , XX( NF ) , X 

FORMAT (IX, 'ERROR IN FINT=R,X OUT OF RANGE, NF , XX ( 1 ) , XX ( NF ) , 

0 X=* ,5X, 13, 3(3X, 010.3) ) 

RETURN 

END 
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Fig. 1: Schematics of a conventional moire system. L.S. is a point light 
source (may be a focused laser beam). M is a parabolic mirror. 6^ 
and G 2 are Ronchi rulings and S is a mat transparent screen. 
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the refracted ray seems to emerge. 
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Fig. 3: Schematic of the experimental setup for postanalyzing phase objects 
G, Ronchi ruling; R.S. and T.S., reference and test signals 
respectively; G.S., gating signal . 
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Fig. 5: Calculated moire intensity vs. phase for various values of f. In each 
figure, the upper two curves are for Aqfa-Gevaert 10E75 plates 
prefogqed with exposure E 0 = 12 erg/cm? and with E'=8 erg/cm2 
and E 1 =18 erg/cm? for the solid line and dashed line, respectively. 

The lower curve is tor Kodak 649F plates with E o =0 and E ' =55 
uJ/cm2. 
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Fig. 7: x - derivative of refraction index for flow simulator - heterodyne 
holographic interferometry. 
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Fig. 8: Object field and the geometry of the transform. 














